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Abstract 

A family of asymptotic solutions at infinity for the system of ordinary 
differential equations is considered. Existence of exact solutions which 
have these asymptotics is proved. 

1 Introduction 

This paper is a reply on the Calogero's question, asked on the NEEDS-2007, 
about justification of an asymptotic solution for the ordinary differential equa- 
tions. The full-dimensional manifold of solutions is the principal interest. Such 
solutions are usually linked with the physically observed phenomena. In the 
case of linear equations an asymptotic expansion of general solution can be con- 
structed by using a fundamental system of asymptotic solutions [51 [15]. We 
present here a similar result for nonlinear equations. 
Main object is the system of differential equations 

— = i fc f(x, t), x G K™, (fc = const > 0) (1) 

considered at t > to > 0. The factor t k is here involved for identification of 
growing right sides, occurring in applied problems; k+ 1 is called Poincare rank. 
Let X(i;a), a G M" be an n-parametric asymptotic solution at infinity. It 
means that the vector function X(i) G C7 1 [io,oo) under substitution into the 
equation fT} yields the residual 

Y(t)^-i fc f(X,i), (2) 
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which is vanishing at infinity: Y(i) = 0(i~ M ), t — » oo, /it = const > 0. We study 
the following problem: is there n-parametric exact solution x(i;a), for which 
the X(i; a) is an asymptotic approximation at infinity: x(t; a) = X(t; a) + 
0{t~ v ), t — > oo, (iv = const > 0)? An affirmative answer is obtained here. 
The proof is reduced to the existence theorem for the remainder. The principal 
feature is the uniformity of the asymptotics with respect to parameter a. 

Justification of asymptotics in general case was not proved up to now and, 
probably, can not be proved at all, although very interesting results are known 
for some types of equations. Results for the linear equations are collected in 
[5J [TS] . Justification of asymptotics for the nonlinear Painleve equations are 
discussed in [21 El E] • Usually existence theorem follows the formal asymptotic 
construction. Various theorems were proved under different conditions which 
depend on the used asymptotic methods [3[ [TJl [Til E2 HU ■ The well known 
Kuznctsov's results [321 [14] deal with the case when X(t) is a formal series 
(of powers, logarithms and exponentials) with constant coefficients. Sometimes 
such series give an asymptotic expansion of general solution [TT] . But often such 
expansion does not correspond to general solution. For example, the pendulum 
equation x + sin x = has few solutions of that type: a single solution x(t) = 
and two one-parametic solutions x(t) = ±ir + 0(exp(— t)), t — ► oo. The other 
solutions are oscillating and belong to two-parametric family. 

Usually the general asymptotic expansions for nonlinear equations are given 
by WKB formulas in the form of power series with oscillating coefficients. 
P.Boutroux studied such two-parametric asymptotics for the Painleve equations 
[2]. In general case the WKB asymptotics depend on n-parameters and just such 
solutions are of interest in physics. For example, existence of two-parametric 
increasing solution of main resonance equations is identified with existence of 
autoresonance phenomenon [H [TU] . Justification of WKB asymptotics for the 
linear equations [5j [16] and for some nonlinear equations 11 can be reduced 
to the Kuznetsov's theorem. But the reduction is unknown and evidently is 
impossible for many equations occurred in physics. Painleves equations studied 
in [H [71 [S] give examples of such type. 

Justification of WKB asymptotics for the Painleve equations is usually based 
on the property of integrability [7] . Nonlinear nonintcgrable equations are stud- 
ied in this paper. We prove an existence theorem and give an estimate for 
the remainder. The presence of the n-parametric solution makes very easy the 
proof. Both the technic of the proof and the derived estimate are determined 
by the properties of the asymptotic solution X(i; a). Similar estimate obtained 
in [5] for the case of single solution X(t) is linked with the properties of matrix 
of linearized system <9xf(X,i). Both of these results are obtained under differ- 
ent sufficient conditions, hence in general case they can not be reduced to each 
other. In order to perform a comparison one have to specify the considered sys- 
tem and the asymptotic solution. We avoid both any asymptotic construction 
and any discussion on the method which yields the asymptotic solution. This 
approach allows us to obtain an existence theorem in general form which may 
be useful for wide range of problems. 

Our estimate is uniform with respect to parameter a and this feature is 
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crucial for the successful application in asymptotic theories. In particular the 
uniformity of the WKB asymptotics proves the orbital stability of the oscil- 
lating solutions. Researches [131 [T41, [H] deal with a single solution which may 
be stable or unstable, see [TT]. As a rule the solution have no any physical 
sense, if it is unstable. Stability of the solution is not discussed in the existence 
theorems [T31 HH H] ■ The problem of stability return us to investigation of n- 
parametric solution which is considered here. Certainly our result does not give 
any information on unstable solutions. 



2 Input conditions 

We consider the equation (JTJ) under assumption that the vector function from 
the right side f(x, t) is smooth (continuously differentiable) in the domain {x £ 
D C M. n ,t > t }. An asymptotics at infinity 

rl d f ( x , *) 

f(x,t),— =0(1), i^oo 

ax 

uniformly on every compact xeKcflis supposed. The principal conditions 
are formulated by means of the asymptotic solution X(i; a). It is assumed that 
the vector function X(£; a) depends on the parameters a = (a\, a n ) and this 
dependence is smooth in some domain a £ Aq C R" . Values of the asymptotic 
solution belong to the compact X(£;a) 6 ICj, C R™ when both t e [t ,oo) 
and parameters a range over any compact: a G A C Aq. Parameters a are 
supposed functionally independent and the Jacobian 



J(t;a) = 




has an asymptotics 

det J(t; a) = t p [a + o(l)], t — *■ oo, (p — const, a — const ^ 0). 

In some sense X(t; a) is a general asymptotic solution. Elements of the Jacobian 
matrix J — {a-ij} have a similar asymptotics (may be different order for various 
number i,j = I, n) 

ai_j(t;a) — 0(t qi ' j ), t — > oo, (q^j — const). 

We shall use the more short formula which may be considered as an asymptotic 
estimate: 

J(t\ a) = 0(t 9 ), t — * oo, (q — maxgij). (3) 
Under these assumptions the inverse matrix has a like asymptotic estimate 

J- 1 (t; a) = 0(f), t^oo. (4) 
lr The case of unbounded asymptotic solution can be reduced to considered one. 
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One can see that the exponents satisfy the relations: p < qn, r < q(n — 1) — p. 

The residual determined by the formula ([2]) depends on the parameters as 
well: Y = Y(t; a). We assume asymptotic estimates: 

BY 

Y(t; a) = 0(t->*), — = 0{t-" +s ), t -> oo, (u,s = const). (5) 
oa 

All asymptotics given above are supposed to be uniform with respect to param- 
eters a G A. 

Remark. All restrictions are induced by examples of WKB asymptotics 
for the Painleve's equations [7J [5] and for the main resonance equations [TU] • 
The exponent /i is a crucial parameter. It corresponds to the length of the 
asymptotic solution. As a rule, /i can be made large as desirecd. While the 
exponents k,p,q,r,s depend on both the considered equations and the family 
of asymptotic solutions. They do not depend on the length of asymptotics /x, 
as is seen in different examples. 



3 Existence theorem 

Theorem 1 Let the vector function X(i;a), a d Aq C R n be an n-parametric 
asymptotic solution of the equation (TTJ) and the properties ^ hold. If the 
residual defined by (J2]) is decreased (O fast sufficiently at infinity: 

p> r + s + 1 and fj, > 2(r + q + 1) + k, (6) 

then for arbitrary compact A C Ao there is a neighborhood of infinity t G [T, oo) 
where the equation (fTJ) has the n-parametric exact solution x(t; a) which has an 
asymptotics 

x(i; a) = X(f; a) + 0(t~ v ), t -> oo, (f = fj,-r-q-l>0) 
uniformly with Va G A. 

Proof 1 The idea of the proof was suggested by M.Fedoruk J^j/ and later was 
applied in [10]. 

We seek an exact solution as the sum x(t;a) = X(£;a) + R(t;a). The de- 
pendence on the parameters a will be omitted for short. The equation for the 
desired remainder R reads 

— =t k [i(X + R,t)-i(X,t)]-Y(t). (7) 

It can be rearranged by identification of the linear part of the operator as follows 

^ = « t [H(t)R + G(R,f)]-Y(t). (8) 



2 In the case of complicated problem the quantity /i depends on the talent of researcher, 
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Here the matrix A4o(t) is calculated by means of derivatives of the known func- 
tions given above: 

ox 

Vector function 

G(R, t) = f(X + R, t) - f(X, t) - <9 x f(X, t)R 

/ias an asymptotics G(R, t) = 0(|R| 2 ) as R — > 0. T/izs asymptotics is uniform 
with respect to t > to- 

M^e desire to demonstrate existence of an exact solution R(i; a) which is 
given in some infinite interval T < t < oo and which is rapidly decreasing at in- 
finity. As usually in such theorems, the differential equation ([8]) mwst oe reduced 
to an integral equation in order to apply the method of successive approxima- 
tions. However, a trivial reduction consisting in integration with respect to t 
does not lead to success because the linear part of the integral operator is not 
contractive in that case. The idea suggested by M.Fedoruk consists in inversion 
of an essential part of the linear operator in the equation jSJ by using the given 
asymptotic solution X(f;a). For this type of inversion we use the Jacobi ma- 
trix 3{t\ a). It is a fundamental asymptotic solution for the linearized equation. 
Indeed, derivation of the starting identity ([2]) with respect to parameters a yields 

^L-t k M (t)J = d a Y. 
dt 

Here the left side represents linearized operator which is applied at J , while the 
right side is the residual which is rapidly vanishing at infinity: d a Y(t; a) — 
0(t s ~L l ), t^oo. 

The Jacobi matrix is used in the substitution R = JC Then the equation 
([8]) is transformed at the equation for the vector function C(t;a): 

^-=M(t)C + g(C,t)-J- 1 Y(t). (9) 
dt 

Main novelty now is the matrix A4. at the linear part 

M(t) = -J~ l d a Y = 0{t r+s ~»), 

which is rapidly decreasing at infinity if the fi is sufficiently large. In particular, 
the norm of the matrix is estimated in W l 

\M(t)\ < M t r+S -", Vt > t , (M = const). 

Nonlinear part is represented by the vector function 

g(c,t) = t k j- 1 G(jc,t). 

It has an asymptotics G(R, t) = 0(|R| 2 ) as |R| — > and this asymptotics is 
differ entiable with R (uniformly with respect to t). Then the inequality 



|G(Ri,t)-G(R 2 ,t)| <Mk(|Ri| + |R 2 |)|Ri-R 2 |, Vt>t 



o 
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take place for all R4, R2 in arbitrary compact IC; here the constant Mr- depends 
on the chosen compact IC. Last relation provides the Lipschitzian property of 
the operator: 

|2(Ci,t) - 0(C 2| i)| < M K t k+r -t«(|Ci| + |C 2 |) • f |Ci - C 2 | (10) 

/or a/Z vector function t 9 d(i)) £ 9 d(i) in every compact IC. 

Differential equation §§§ is equivalent to the integral equation 

C(i)- f\M( V )C(r 1 ) + g(C(7 1 ),7 1 )}dn = Z(t). (11) 
TTie vector function from the right side 

/oo 
J- l Y{rj)d n 

is well known and has an asymptotics Z(t) = 0(t r+ i — » 00. We are looking 
for a solution C(t) which has a like asymptotic behavior. 

In order to prove the theorem, we have to introduce some Banach space. An 
appropriate space of continuous vector functions C(t) € C\[T, 00) is determined 
by the weight norm 

||C|| A = sup* A |C(*)|, X = n-r-l. 

t>T 

The boundary T — const > to will be chosen from requirement of contractiveness 
of the integral operator in what follows. 

The first condition from (jJJJ) taken in the form A > q and the inequality (flQ|) 
provide the estimate 

\g(d,t) - G(c 2 ,t)\ < M K t k + r + 2 *- 2 * ■ (Hdiu + ||c 2 || A )||c a - c 2 || A (12) 

which is valid for all vector functions Ci(i),C2(i) in arbitrary compact IC C 
C\[T, 00). If we take here integral with t, then we can derive the Lipschitzian 
inequality 



f mC^-G^v)] dr, 

J OO 



^ < M K 2KT k+r+2q+1 - x ■ ||d - d|U = L K -\\C X - d|U, 



(13) 

which is valid for all ||Ci||a, || C2 1| A < K- The Lipschitzian constant 

L K = M K 2KT k+2 ^ +q+1 ^-^ 

can be chosen arbitrary small (for example, Lk < 1/2) under large T, because 
the exponent k + 2(r + q + 1) — (i < is negative because of the second condition 
(|6|). The relation (|13|) provides contractiveness of the nonlinear part of the 
integral operator. 
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The norm of the linear part of the integral operator can be estimated by 
integration of the matrix norm M. (t) : 



supi A / n~ x \M{n)\dn < M x supt x 



= M 1 T r+s+1 -' 1 /(r + s - n - X + 1) = Lq. 

The constant Lq can be made small as desired (say, L < 1/2) under assumed 
condition fi > r + s + 1. It is enough to choose an appropriate large T . The Lq 
is used as a Lipschitzian constant for the linear part of integral operator: 



M(r,)[C 1 (v)-C 1 (r 1 )}dr 1 



x <L \\C 1 -C 2 \\ x . 



Thus, if the value of T is sufficiently large then the integral operator from 
the equation turn out to be contractive in Banach space C\\T, oo). The 

right side of the equation Z(f) = 0(i r+1 ~ M ), t — > oo belongs the same space 
Z(t) £ C\[T, oo) as X = /jt — r — 1. Then the fixed-point theorem can be applied to 
the integral equation (jlip . It implies existence of the solution C(t) S C\[T, oo). 
This vector function is differ entiable because of the identity As the equation 
© is equivalent to 0, then R(t) = JC(t) is a desired solution of the residual 
equation ([7]) and it has an asymptotics 

R(f) = JC(t) = 0(t q - x ), t^oo, (A - q = v). 

It must be point out that all found estimates are uniform with respect to 
parameters a in arbitrary compact A. So the proof, as given above, is suitable 
for the both case either fixed parameter a or fixed compact: a E A C Aq . In the 
last case the estimate of the remainder 

H(t; a) = JC(t; a) = 0{r v ), t -> oo 

is uniform with respect to V a G A. This completes the proof. 

Remark. It may be occurred in application at specific problems that the 
order of the asymptotic estimate of the remainder R(t) = 0(t^~ r ~ q ~ 1 ) is very 
rough. It is lack of the theorem given for general case. More precise estimations 
can be derived, if a full asymptotic solution in the form of infinite asymptotic 
series is known; see |10j . 
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